Scalar interaction limits from the (3-v correlation of trapped radioactive atoms 
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We have set limits on contributions of scalar interactions to nuclear /3 decay. A magneto-optical 
trap (MOT) provides a localized source of atoms suspended in space, so the low-energy recoiling 
nuclei can freely escape and be detected in coincidence with the (5. This allows reconstruction of 
the neutrino momentum, and the measurement of the (3-v correlation, in a more direct fashion than 
previously possible. The f3-i> correlation parameter of the + — ► + pure Fermi decay of 38m K is 
a=0.9981±0.0030tp;oo37, consistent with the Standard Model prediction 5=1. 
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The angular correlation of neutrinos and beta's in nu- 
clear beta decay is historically one of the main experi- 
mental probes of the vector and axial vector nature of 
the weak interaction 0. A recent experiment using 32 Ar 
decay is the only (3-v correlation measurement in pure 
Fermi decay, which is sensitive to scalar interactions 0]. 

We use a magneto-optical trap (MOT) Q to provide 
a backing-free source of atoms with well-defined position 
and negligible thermal energy. We then detect the low- 
energy nuclear recoil in coincidence with the emitted (3 + , 
and directly deduce the v direction and the (3 + -v corre- 
lation. We also determine critical response functions of 
our detectors in situ from the decays themselves. Atom 
trap (3-v experiments are also pursued elsewhere Q. 

In the + — > + Fermi decays the leptons carry away 
no net angular momentum. Back-to-back (3-v emission is 
forbidden in the Standard Model, because the W vector 
boson exchange produces leptons with opposite helicity 
and their spins add to one. The angular distribution is 



W(9 Pv ) = l + b^ + 



a— cos(9/3v). 



The (3-v coefficient a is +1 for W exchange, and a is -1 for 
a scalar boson producing same-helicity leptons. In terms 
of scalar coupling constants Cs and C' s , and assuming 
=1 !: 



for simplicity Cv=C' v - 
a=[2-(\C s \ 2 



\C's\ 2 )}/(2- 



|Cs| 2 + |^| 2 ), 



b = -2Vl - (aZ) 2 Re(C s + C' s )/{2 + \C S 



\C'J 2 ) 



The limit on the Fierz interference term b from the depen- 
dence of + ^0 + decay strengths on (£73) is very strin- 
gent, 6=-0.0027±0.0029 @], but the coupling C s +C' s 
describes scalars that couple only to the left-handed v. 
Measurements of a constrain scalar interactions indepen- 
dent of chirality or time- reversal properties 0, 0] • 



The isobaric analog decays of the pure Fermi tran- 
sitions are well characterized. Lowest recoil-order cor- 
rections to the allowed approximation value of a=l 
(<3xl0~ 4 in our case, 38m K) do not depend on nuclear 
structure, and higher order corrections are <0.0002 
Radiative corrections (see below) can also be calcu- 
lated to the order required independent of nuclear struc- 
ture. In addition, 38m K decay is known to proceed 
cleanly to the ground state, with experimental limits on 
excited-state branches of < 2x10 Disagreement 
with a=l greater than these corrections would be from a 
standard model-violating scalar interaction. 

A scalar term could be produced by the exchange of 
scalar bosons found in standard model extensions Q. 
There is a phenomenological window open for sleptons, 
although couplings to the first generation are thought 
to be small p}. A QCD-induced scalar interaction is 
a second-class current that is <5xlCP 5 in the Stan- 
dard Model /3-v correlations are complementary to 
scalar constraints from pseudoscalar n^ve decay |l2j . 
The possibility of constraining scalar interactions from 
loop corrections to v masses is considered in Ref. [l3j . 
The best previous (3-v experiment used /3-delayed pro- 
tons from 32 Ar decay to determine a=a/ (1 + bmp/(Ep)) 
to be 0.9989±0.0052±0.0036 @. 

Our apparatus was also used to constrain massive v x - 
v e admixtures Q. A beam of 38m K ions (t 1 / 2 = 0.924 s, 
= 5.022 MeV) is produced at TRIUMF's ISAC fa- 
cilities dSj, sto ppe d and released as neutral atoms with 
a 900°C Zr foil and captured with w 10~ 3 efficiency 
in a vapor-cell MOT. The MOT traps only the 38m K 
and none of the ground state 38 K beam contaminant . To 
escape backgrounds from untrapped atoms of 38 K and 
38m K, we transfer the trapped atoms with 75% efficiency 
by a chopped laser push beam to a 2 nd MOT equipped 
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FIG. 1: Top view of the 2 nd MOT apparatus with the recoil 
and (3 detectors. 

with the nuclear detectors (Fig. QJ. The duty cycle en- 
tails: push atoms from the first trap for 20 ms; wait 50 
ms to transfer; change the 2 nd MOT laser frequency and 
power to minimize cloud size; wait 1 ms to let cloud reach 
equilibrium; count for 150 ms; repeat [r7|. No atoms are 
lost from the trap during the frequency switch. The MOT 
force uses laser light and a weak (dB z /dz = 20 G/cm) 
magnetic quadrupole field, so the Ar recoils escape the 
trap without perturbation. We accelerate the positive 
Ar ions produced by electron shakeoff 0] with a uni- 
form electrostatic field to separate them in time-of-flight 
(TOF) from the neutral Ar° atoms. 

The /3-telescope is a position-sensitive 22x22x0.49 mm 
double-sided Si-strip detector (DSSSD) backed by a 
06.5 cm x 5.5 cm long BC408 plastic scintillator. The 
telescope is separated from the trap vacuum by a 125^m 
thick Be foil located 2mm from the DSSSD to minimize 
angle straggling. The telescope coincidence rejects 95% 
of the 2.17 MeV 7-rays from untrapped 38 K ground state. 
The gain is actively stabilized at the low count rates of 
< 200 Hz using a stabilized light pulser. 

The Ar recoils, which have 0-430 eV of initial ki- 
netic energy, are detected by a Z-stack of three uncoated 
microchannel plates (MCP). A fixed aperture defines a 
24.0 mm active diameter for the TOFfE^] analysis (see 
below). The resistive anode position readout is calibrated 
with a mask and an a-source to have 0.25 mm resolution 
within the 20 mm diameter used for the reconstructed an- 
gular distribution analysis. The E field accelerates the 
Ar +1 ions to 4.8-5.3 keV. We measured the MCP effi- 
ciency in this energy range to be constant to accuracy 
0.0060 by comparing the rate of /3-recoil coincidences for 
four values of E$. The (3-v correlation analysis is done 
with the ions, because the efficiency for neutral recoils 
is not as well understood. The angle dependence of the 
MCP efficiency was assumed constant over the small im- 
pact angles of ±5°, with error (Table HJ) spanning the 



small effect seen in the literature [19J , consistent with our 
analysis of recoils that uniformly illuminate the MCP. 

We maintain a population of w 2,000 atoms of 38m K in 
the detection MOT. The trap lifetime, limited by resid- 
ual gas, is 45 s, so 97% of the 38m K atoms decay while in 
the trap. Atoms on the walls produce a j3 + singles back- 
ground of < 2% and a negligible coincidence background 
consistent with accidental coincidences, measured by de- 
liberately releasing the atoms to the walls. Ions from the 
walls are excluded from the MCP by the electric field. 
Ions from the trap strike no material before reaching the 
MCP. The E field electrodes are made from glassy carbon 
to minimize [3 + scattering effects. 

The average trap-MCP distance is determined to be 
61.08 ± 0.01 mm from a fit to the leading edge of the 
TOF peak of the fastest Ar° recoils (Fig. EJ)- These were 
shown to be Ar° by applying E fields of 400 and 800 
V/cm; the leading edge was undistorted by any detection 
of the r e = 260 ns |2fJ Ar~ metastable state. 

We image the cloud by photoionizing a small fraction 
of the 38m K atoms with a pulsed laser and accelerating 
them to the MCP. The x, y, and z distributions (see Fig. 
1) are fit well with Gaussians of 0.8, 1.1, and 0.65 mm 
FWHM. The z distribution (along the trap-MCP axis) 
limits the timing resolution for Ar +1 recoils to 5 ns. Two 
CCD cameras image the trap laser fluorescence, and the 
trap centroid was kept constant to ±0.05 mm. 

We performed two independent analyses of the data 
set. In the first analysis, we fit the TOF spectra of ion re- 
coils for various (3 + energy cuts (Fig. |3J to a Monte Carlo 
(MC) simulation based on GEANT |2l). Qualitatively, 
for fixed E^, the recoil TOF increases monotonically with 
cos^,,). The TOF of the ions with most sensitivity to 
a increases with decreasing Ep. We fit simultaneously to 
A. r +i,+2,+3 charge states. 

In the second analysis we use the complete momentum 
information measured for the [3 and the recoil to deduce 
the momentum of the v and the f3 — v angle (Fig. [5J. The 
kinematics are over-determined for p rC coii< Qp+/c |22j, 
so this was done either using the measured Ep or de- 
termining it from the recoil momentum. The measured 
angular distribution is fit to the MC simulation as a func- 
tion of a, and agrees with the TOF[E^] analysis. 

We present the details of the TOFfE^] analysis and 
detailed evaluation of systematic errors. This analysis 
lets us constrain critical physical and instrumental ef- 
fects, but requires excellent j3 telescope characterization. 
Table H] shows systematic errors, determined by MC sim- 
ulations varying each parameter by its possible error and 
determining its effect on a, with other parameters refit 
as appropriate. Some errors in the table are summaries 
of more than one correlated systematic error. These line 
items are uncorrelated, so we add them in quadrature to 
determine the total systematic error. 

We can test the MC simulation of the (3 + energy line- 
shape with the /3-recoil coincidences. From Fig. [2J the 
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FIG. 2: (color online) Bottom: Scatter plot of recoil TOF 
vs. Tfj with one dot shown for each of 500 000 events. The 
suppressed back-to-back lepton emission produces longer re- 
coil TOF. The E field separates the Ar charge states. The 
analysis cuts are shown. Top: TOF projections of the 2D 
scatter plot. The 0.1% background at TOF«1020 ns is from 
/3's scattering off the MCP into the /3 telescope and can be 
rejected kinematically. 

(E) field/trap width: 0.0017 
E field nonuniformity 0.0014 
~Eif3+ Detector Response: 

Lineshape tail/total 0.101T0.006 0.0006 
511keV Compton summing/total to 10% 0.0009 
Calibration including nonlinearity 0.0017 
MCP Eff[E A r+] measured constant 4.8-5.3keV 0.0007 
MCP Eff[0]/XY trap position 0.0008 
e~ shakeoff dependence on p rC coii s = 07_° 014 to.oois 

TABLE I: List of a uncorrelated systematic errors. 



Ep spectra in coincidence with neutral recoils for TOF 
intervals from 1500 to 1800 ns are peaks determined by 
the detector resolution and the angular acceptance, and a 
tail determined by /3 + backscatter, bremsstrahlung, and 
^20% of the tail from scattering off inactive volumes. 
From this and from detailed kinematic reconstruction of 
the /3-Ar +1 coincidences, we have determined that both 
the size of this tail and the 511 keV Compton summing 
agree with the MC simulation; they are listed in Table I. 

The Ep calibration is determined by a MC fit to 
the energy spectrum in coincidence with recoils with 
370 ns < TOF < 900 ns, which includes all the 
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FIG. 3: (color online) Top: lowest of 16 Tp bins for the MC 
fit to the Ar +1,+2,+3 data, and residuals. The confidence level 
for the entire fit is 52%. Data has been binned to show sensi- 
tivity to a. Fits are done with 4 ns bins. The dip in Ar +1 is 
from the finite MCP size: the dashed curve has an artificially 
larger MCP collecting all ions. Bottom: Fit to reconstructed 
angular distribution, showing residuals. Also shown are the 
effect of a change in a by 0.005, and the effect of recoil energy- 
dependent shakeoff. 



observed Ar +1,+2,+3 recoils. We use the expression 
xadc — xq + C2 Tp/(l+qTp), with nonlinear term 
q = (0.33±1.49)xl0~ 3 MeV -1 . The calibration param- 
eters are not sensitive to the value of a. Use of an E_a 
calibration determined from the j3 singles spectrum over 
the fit range produces a value of a consistent within the 
error in Table I. The fit to the coincidence energy spec- 
trum has x 2 /N=21.8/23, and the fit to the singles en- 
ergy spectrum has x 2 /N=10.8/ll. We use the Fermi 
function and corrections of |2^|. 

We are working to extend the experiment to Tp <2.5 
MeV to independently determine b. The fits become 
poorer, and systematic errors from (3 scattering, the low- 
energy lineshape tail, and possible additional sources are 
not fully understood. The Tp cutoff eliminates almost all 
backscattered events, as well as all possible contamina- 
tion from untrapped 38 K ground state decays. 

We include order-a radiative corrections. These are 
dominated by undetected momentum carried away by 
real bremsstrahlung photons, which we include in the 
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MC event generator|24| . They change a by 0.003 in the 
32 Ar experiment Because the /3 + energy spectrum 

is also affected, and we use the (3 + spectrum itself for our 
Efj calibration, our net result is that the radiative cor- 
rections change a by considerably less in our experiment. 

Three independent measures determine the E field. 
The leading edge TOF of the Ar ion spectra implies 
-Es=807.7±0.16 V/cm, independent of a and consistent 
for all charge states. The field nonuniformity is con- 
strained by the TOF of the photoionized 38m K atoms, 
and by a population of 'wrong-way' recoils produced from 
f3's firing the MCP, which give central values 807.7 and 
808.3 V/cm. The nonuniformity is < 1.0 V/cm/cm and 
the resulting error in a is 0.0014. 

We collect 89%, 99.6%, and 99.9% of the Ar +1 <+ 2 <+ 3 
ions in coincidence for Tp > 2.58 MeV. This finite accep- 
tance is a source of systematic error (see Fig. 3); since a 
fixed aperture defines the MCP size, the acceptance con- 
tributes part of the dependence of a on the E field and 
trap position, as quantified by the MC analysis (Table I). 

Dependence of the probability of electron shakeoff on 
the recoil ion energy has been seen in 6 He /3~ decay [lif . 
A recent simple estimate relates this effect to oscillator 
strengths and suggests that it is larger in (3 + decay 
The recoil energy spectrum to lowest order is distorted 
by (l+sE rcc ). We constrain this effect experimentally 
by fitting s and a simultaneously in our TOFpSg] fit 
for Ar +1,+2,+3 . We only include s in the Ar +1 spec- 
trum, because the model of Ref. [2{| using semiempir- 
ical oscillator strengths [i^ suggests that s for Ar+ 2 (or 
Ar+ 3 ) would be 0.11 (or 0.05) the size of s for Ar +1 . 
We find s=— 0.013±0.020, a result in a nonphysical re- 
gion with one a upper limit s < 0.014 and change in a: 
A a = 0±q 0018 . The estimate of Ref. JU is s=0.031. A 
similar fit of the Ar +1 reconstructed angular distribution 
to a and s gives consistent result and error (Fig. . We 
can constrain s and a simultaneously because the greatest 
sensitivity to a is at the null in the angular distribution, 
and because we fit as a function of E^. A fit to the to- 
tal TOF spectrum summed over all E^ would be more 
strongly correlated with the recoil momentum spectrum. 

Our fit values for a and b are strongly correlated in 
the Ep region used. Although we fit as a function of Ep, 
careful investigation of the correlations shows that the 
physical observable we report here is, for |6|<0.04, effec- 
tively indistinguishable from that reported by Ref. Q, 
a= a/(l + bmp/(Ep)), but with (Ep)=3.3 MeV. We find 
5 = 0.9981 ± 0.0030±°-°037> in agreement with the Stan- 
dard Model. If we vary b from -0.0075 to +0.0021, the 



90% confidence range of the Fierz interference term limits 
in Ref. [g, then a changes from 0.9971 to 0.9984, while a 
changes by <lxl0 -4 . Our measurement has comparable 
errors to Ref. 0] with an entirely different experimental 
method. 

We acknowledge TRIUMF/ISAC staff and comments 
by A.R. Young, R.M. Woloshyn, J. Ng, B.A. Campbell, 
LA. Towner, and P. Herczeg. Supported by the Na- 
tional Research Council of Canada through TRIUMF, 
the Natural Sciences and Engineering Research Council 
of Canada, the Israel Science Foundation, and WestGrid. 



* Deceased 

[1] E.D. Commins and P.H. Bucksbaum, Weak Interactions 
of Leptons and Quarks, Cambridge University Press, 
1983, Fig. 5.3. 

[2] E. G. Adelberger et al, Phys. Rev. Lett. 83 1299 (1999); 
erratum 83 3101 (1999). a is being reevaluated: see A. 
Garcia, Nucl. Phys. A746 (2004) 298c. 
[3] E.L. Raab et al, Phys. Rev. Lett. 59 2631 (1987). 
[4] N.D. Scielzo et al, Phys. Rev. Lett. 93 (2004) 102501. 
[5] J.D. Jackson, S. B. Treiman and H.W. Wyld, Phys. Rev. 

106 (1957) 517; Nucl. Phys. 4 (1957) 206. 
[6] I.S. Towner and J.C. Hardy, J. Phys. G 29 (2003) 197. 
[7] P.Herczeg, Prog. Part. Nucl. Phys. 46/2 413 (2001). 
[8] B. Holstein, Rev. Mod. Phys. 46 789 (1974). 
[9] F. Gliick, Nucl. Phys. A628 493 (1998). 
[10] E. Hagberg et al. Phys. Rev. Lett. 73 396 (1994) 
[11] B.R. Holstein, Phys. Rev. C 29 623 (1984). 
[12] B.A. Campbell and D.W. Maybury, Nucl. Phys. B 709 

419 (2005). 

[13] T.M. Ito and G. Prezeau, hep-ph/0410254 (2004). 
[14] M. Trinczek et al, Phys. Rev. Lett 90 012501 (2003) 
[15] M. Dombsky et al, Rev. Sci. Instrum. 71, 978 (2000). 
[16] D. Melconian et al, Nucl. Instru. Meth. A 538 (2005) 
93. 

[17] T.B. Swanson et al, J. Opt. Soc. Am. B 15, 2641 (1998). 

[18] T.A. Carlson, Frances Pleasonton, and C.H. Johnson, 
Phys. Rev. 129 (1963) 2220. 

[19] G.W. Fraser, Int. Jour. Mass Spec. 215 13 (2002). 

[20] I. Ben-Itzhak et al, Phys. Rev. A 38, 4870 (1988). 

[21] GEANT version 3.12, CERN (1994). It includes the low- 
energy physics of EGS4. 

[22] O. Kofoed-Hansen, Dan. Mat. Fys. Medd. 28 no. 9, p. 1 
(1954). 

[23] D. Wilkinson, Nucl. Instru. Meth. A290 509 (1990). 
[24] F. Gliick, Comp. Phys. Comm. 101 (1997) 223. 
[25] N.D. Scielzo et al, Phys. Rev. A 68 (2003) 022716. 
[26] D.A. Verner et al, Astrophys. J. 465 487 (1996). 



